Microscopic theory of Cooper pair beam splitters based on carbon nanotubes 
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We analyze microscopically a Cooper pair splitting device in which a central superconducting 
lead is connected to two weakly coupled normal leads through a carbon nanotube. We determine 
the splitting efficiency at resonance in terms of geometrical and material parameters, including the 
effect of spin-orbit scattering. While the efficiency in the linear regime is limited to 50% and decay 
exponentially as a function of the width of the superconducting region we show that it can rise up 
to ~ 100% in the non-linear regime for certain regions of the stability diagram. 



PACS numbers: 73.63.-b, 74.45. -|-c, 73.63.Fg 

Introduction: Producing entangled electron pairs in a 
solid state device from the splitting of Cooper pairs [T] is 
a challenging possibility which is starting to generate an 
intense experimental effort [5H1] . A basic splitting device 
is a three terminal system with a central superconduct- 
ing lead (S) in between two normal (N) ones as depicted 
in the upper panel of Fig. [l] When a Cooper pair is 
injected from the S lead it can either be transmitted as a 
whole to one of the N leads by means of a local Andreev 
process (Fig. [I^) or split so that each of the electrons in 
the pair is transmitted to a different lead, which corre- 
sponds to a cross erf Andreev process (CAR) (Fig. ^) [5]. 
While initial experimental devices were based on nano- 
litographically defined diffusive samples more recent 
experiments are oriented towards tunable double quan- 
tum dots systems based either on carbon nanotubes [3] 
or InAs nanowires [4]. In spite of the difference in mate- 
rials the systems realized in both experiments were con- 
ceptually equivalent. They did correspond, however, to 
different physical regimes: while in Ref [5] the hybridiza- 
tion by direct tunneling between the dots was dominant, 
in Ref. [4] the direct tunneling appeared to be negligi- 
ble. In both works the Cooper pair splitting action was 
demonstrated indirectly by analyzing the changes in the 
behavior of the conductance when going from the normal 
to the superconducting state. Both works pointed out an 
unexpectedly high efficiency for CAR, much higher than 
what would be predicted by theories which do not take 
into account the direct inter-dot tunneling [T]. Still fur- 
ther experimental and theoretical efforts are needed in 
order to demonstrate the splitting unambiguously and to 
reach the nearly 100% efficiency which could be necessary 
for entanglement detection ^6j. 

In the present work we analyze microscopically the case 
of double quantum dots (DQD) defined on single-walled 
carbon nanotubes (SWCNTs) and show that the two 
regimes of Refs. H] can be reached in metallic or semi- 
conducting tubes. We consider the situation illustrated 
in Fig. [T]; and[TJi where the central electrode modifies the 
electrostatic potential and induces a pairing amplitude on 




FIG. 1: (Color online) Schematic representation of local (a) 
and non-local (b) Andreev processes in a generic Cooper pair 
splitter, (c) Specific geometry considered in this work: a finite 
SWCNT coupled to normal leads at its ends and a central 
superconducting lead. The lower image (d) illustrates the 
potential profile along the tube. 



the portion of the tube underneath without breaking its 
continuity. In agreement with the experimental obser- 
vations we show that in this case the splitting efhciency 
decays rather weakly with the width of the central elec- 
trode [1]. Our results also suggest how to increase the 
splitting efficiency up to a level close to 100% by operat- 
ing the devices in the non-linear regime. We begin with 
a SWCNT in the normal state without e-e interactions in 
order to analyze the inter-dot coupling. Subsequently, we 
switch on superconductivity in the central electrode and 
study the probability of the CAR processes and the split- 
ting efficiency. Finally, we use a minimal model to an- 
alyze the effect of the electron-electron interactions and 
the non-linear regime, where the efficiency can be highly 
enhanced. 

Basic modeling: We focus on zig-zag SWCNTs which 
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allow to consider both metallic and semiconducting cases 
[8]. If the coupling to the central lead is sufficiently 
smooth on the atomic scale we may assume that inter- 
valley scattering is weak and the K-K' degeneracy is pre- 
served. For this case and when the radius is of the order 
of 1 nm or smaller it is important to consider curvature 
effects which produce a finite band gap in metallic tubes 
and enhances the effect of spin-orbit (SO) interactions. 
We use two complementary approaches for describing the 
electronic states in the zig-zag SWCNT: a tight-binding 
(TB) model and a continuous description based on the 
Bogoliubov-de Genes-Dirac equations. While the latter 
allows analytical insight, the numerical TB calculations 
allow to identify effects due to deviations from linear dis- 
persion, disorder or arbitrary spatial variation of the elec- 
trostatic potential along the tube. 

Within the continuous description the system is char- 
acterized by the equations 
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where H!^ ^ — —ihvpdx ■ ty^ + Thvpqnay + tSos — rSisay + 
V{x) is the normal state effective Hamiltonian for the 
n mode (corresponding to a quantized momenta 
around the tube), A{x) is the induced pairing ampli- 
tude and V{x) the electrostatic potential profile along 
the tube. In these equations are Pauli matrices in 
sublattice space, and r, s = ± correspond to the val- 
ley and spin indexes respectively. Finally, the terms in 
Sq and 6i take into account the SO interaction as in 
Refs. [inilTT]. The quantized momenta take the values 
= 1) ~ Icurv, with p = A^modS = 0, ±1 

and N being the number of atoms in the cross sec- 
tion. Depending on whether p = Oorp = ±1 the 
tube is metallic or semiconductor, respectively. The cur- 
vature effect is included in Qcurv = E^urv l^vp, where 
Ecurv ^ Tr'^\Vpp^\/4N'^ with \Vpp^\ ~ 2.7eV. In the nor- 
mal homogeneous case the corresponding energy levels 
for longitudinal wavevector k are thus given by EJ^ g{k) — 

1wf\/ {qn + Tsdi)^ + k'^ +SoTs. The transport properties 
can be expressed in terms of Green functions which sat- 
isfy {E — 'Hr,six))Gr,s{x, x') — 5{x — x'), whcrc Hr.six) 
denotes the full Hamiltonian on the left hand side of Eq. 
([l]). We obtain these quantities by solving first for the 
uniform finite regions and then matching the result using 
the method of Ref. [7^ . 

Normal state: We start by analyzing the linear con- 
ductance along the tube, Glu, when the central lead is 
in the normal state. In Fig. [2] we show a map of Glb. in 
the VgL — VgR plane, obtained using the TB model in the 
usual nearest neighbors approximation with a hopping 
parameter t = Vppjr ■ As a first approximation the poten- 
tial profile along the tube is assumed to change discontin- 
uously as represented by the dashed lines in Fig. [TJi. The 
lateral leads are modeled by ideal ID channels weakly 
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FIG. 2: (Color online) Left panel: Conductance map for a 
metallic tube (A'^ = 12) in the normal state for the p-n- 
p region with (lower) and without (upper) SO interactions. 
Right panel: Same for the semiconducting tube (A'^ — 11) 
but in a logarithmic scale. The geometrical parameters are 
W = Wl,r = nOnm. 



coupled to the ends of the tube, as represented schemati- 
cally in Fig. [ij;. We fix the tunneling rates to these leads 
to a value ^ O.Oli which is consistent with the con- 
ductance values observed in Ref. [HI for the lowest energy 
states of a SWCNT quantum dot. To model the effect of 
the central lead we rely on the results of ab-initio calcu- 
lations for the case of Al electrodes [TH [T^ . According 
to Ref. [I2] these produce a n-doping effect, leading to 
a shift of the tube bands Eps ~ —0.5eV for an ideal in- 
terface. On the other hand in the normal state it would 
induce a broadening of the tube levels of the order of a few 
meV [13] which suggests a typical value Fg ~ ImeV for 
the corresponding tunneling rate. As in the experiments 
of Ref. [3] the length of the central region is set initially 
to ^ 200nm. We consider tubes with N = 11, 12 which 
corresponds to radii R ~ 0.43, 0.47nm for the semicon- 
ducting and metallic cases respectively. In the metallic 
case curvature effects lead to the opening of a narrow 
gap, which can be estimated as Eg ~ Ecurv — -ibmeV. 
The curvature gap is apparent in the upper left panel of 
Fig. [2] On the other hand, in the semiconducting case 
the gap is Eg ~ 412'meV (top right panel of Fig. [2]). It 
should be noticed that for these diameters and for gate 
potentials of the order of 0.5eV only the lowest energy 
mode, corresponding to rt = 0, gives a significant contri- 
bution to the transport properties of the tubes that we 
discuss below. 

For positive VgL, VgR, i.e. in the p-n-p regime the con- 
ductance displays a DQD behavior as shown in Fig. [2] 
The metallic case (left panels) exhibits an anticrossing 
pattern similar to the one found in the experiments of 
Ref. ^ . As the gate potentials VgL , VgR become increas- 
ingly positive the conductance map exhibits resonances 
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along lines VgL + Vgn ~ const indicating the delocaliza- 
tion of the electronic states due to Klein tunneling. The 
confinement of the dot states is much more pronounced 
in the semiconducting case where the Klein tunneling is 
less significant. We have used a logarithmic scale in this 
case in order to enhance the visibility of the conductance 
peaks. 

When SO scattering is introduced (lower panels in Fig. 
[2]) there is a general splitting of the conductance peaks 
of the order of ~ 2meV due to the breaking of the spin- 
valley degeneracy. Close to the gap edges this splitting 
is of the same order as the mean level separation. 

Superconducting state: When superconductivity in the 
central lead is switched-on pairing correlations within 
the tube are induced by proximity effect. The size of 
the induced gap is set by (i.e. of the order of 
ImeV). We shall assume that temperature is zero and 
that the energy E of the injected electrons from the nor- 
mal leads is smaller than A^. Then Ral{E) and Rar{E) 
denote the local Andreev reflection probabilities at the 
L,i? leads while Tcar{E) correspond to the CAR pro- 
cesses. When operated as a beam splitter a finite voltage 
difference V is applied between the S and the N leads 
and the non-linear conductance is given by Gl(j?)(1^) = 

G(,{Tcar{V) + TcARi-V) + 2RALiR}iV)), with Go = 

2e^//i (notice that at finite energy in general Tcar{E) ^ 
Tcar{—E) due to the breaking of the electron-hole sym- 
metry). Thus we can define the splitting efficiency as 
77 = Go[Tcar{V) + Tcar{-V))/{Gl{V) + Gr{V)). 

Within both the Dirac and TB models it is found that 
the CAR coefficients decay exponentially on the scale 
5(g) = hvp/K^Jl - {hvfq/EpsY, where q = qQ±Si/hvp 
exhibiting oscillations on the scale Xp — fivp /\Eps\, as 
illustrated for the linear regime, V = 0, in Fig. |3] In 
these plots we have fixed the gate voltages at the values 
indicated by the circles in Fig. [2] The CAR probability 
decays more slowly in the metallic case {N = 12) due to 
the longer effective coherence length ^(g) of this case. In 
both cases, however, the decay is remarkably slower than 
in the case of a 3D bulk BCS superconductor where the 
prediction is ~ exp{—2W/^)/{kpW)'^ [H], indicated by 
the dashed lines in Fig. |3] This is a consequence of the 
single channel character of the connection between the 
dots in the present system and explain the rather large 
efficiency values estimated in recent experiments As 
can be observed the efficiency i] decreases from 0.4 at 
W ^ 200nm to < 0.05 at ~ 70Qnm in the semicon- 
ducting case, while it varies between 0.5 and 0.2 for the 
metallic tube within the same W range. For sufficiently 
large W the overall evolution of ij is well described by 
the expression 77 ^ l/(l + exp2VF/5(g)). This qualitative 
behavior is also found for a smoother potential profile 

m- 

Coulomb interactions and non-linear regime: For an- 
alyzing the effect of electron-electron interactions we 
first map the system into a minimal model in which we 
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FIG. 3: (Color online) Evolution of the CAR probability 
(blue) and the splitting efficiency rj (red) in the linear regime 
as a function of the length W of the central electrode for a 
metallic SWCNT with A*' = 12 (upper panel) and a semicon- 
ducting one with A = 11 (lower panel). The gate potentials 
VgL,gR are fixed at the points indicated by the circles in Fig. 
[2] The dashed lines represent the decay of the CAR prob- 
ability for a 3D bulk superconductor multiplied by a factor 
1000. 



keep just one two-fold degenerate electron level -El,_r in 
each dot [IQ. In the combined dot-Nambu space the 
model properties can be expressed in terms of bispinor 
fields = (dp^, dj^^) where ^ ^ L,R and c?j,g. creates 
dot electrons. In the absence of interactions this re- 
duced model is characterized by a retarded Green func- 



tion matrix of the form = 



E-ha+it-t{E) 



where (ft.o)^!',^^ = E^5f^^5api-l)°'~^'^ , with ^i,v = L,R, 
a,/3 = 1, 2 arc the Nambu indexes, (F)^(^.q^ 



with T^ = Tfj^ao/W^, correspond to the effective tunneling 
rates to the normal leads and S is a matrix self-energy 
describing the coupling with the central superconducting 
region. The whole analytical expression for E is given in 
the supplementary information. 

Interactions are introduced by assuming a constant 
charging energy Ul,r ^ A^ acting on each dot level. 
We take them into account within the equation of mo- 
tion (EOM) technique with a Hartree-Fock decoupling 
at the level of the two-body Green functions [T3] . This 
approximation is valid when Kondo and exchange corre- 
lations between the dots can be neglected. Further sim- 
plification is achieved in the limit [/^ ^ 00 where we 
I -1 



find G 



-i+iF-E 



with .g= (S-/io)~Ml-^c 



and (Aao) 
mean values n,. 



- n^S^iv^ajS- The evaluation of the 



must be performed self- 
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consistently. The relevant transport coefficients are fi- 
nally computed as Tcar{E) ^ ^rLi^R\GLRME)\'^ and 

RaL{R}=4Fl(B.)\^LL{RR).12\'^- 
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FIG. 4: (Color online) Color map of the splitting efficiency 
within the minimal model with parameters corresponding to 
a semiconducting tube with W ~ 700nm (indicated by the 
arrow in Fig. [sf in the linear (a) and non-linear (b) regimes. 
In the latter case the dot levels are varied along the line 
El ~ —Er indicated by the dashed red line of (a). The 
white dashed lines indicate the maxima in the spectral den- 
sity which is shown in (c) for the two dots along this line. 
Local and non-local Andreev processes at finite V are indi- 
cated by the black arrows. 

The main effect of interactions within this approxi- 
mation is to shift the resonances and to reduce their 
width, roughly as (1 — n^)T^. Then, the CAR and 
the local Andreev probabilities are reduced by a factor 
(1— ni)(l— ^2) and (1 — n^)^ respectively, which therefore 
does not modify significantly the efficiency at resonance. 
The color map in Fig. |4|i shows the efficiency in the 
linear regime corresponding to the semiconducting case 
with W ~ 700nm (arrow in Fig. [s]) and for the region of 
gate voltages indicated by the circle in the right panel of 
Fig. [2| As can be observed, rj exhibits maximum values 
at the crossing point between the resonances of the order 
of 0.1 which are slightly higher than the values found in 
the non-interacting case. The efficiency reaches a maxi- 
mum of the same magnitude along the line El ^ —Ej^ 
(red dashed line) . What is much more remarkable is that 
the efficiency along this line can rise up to 100% in the 
non-linear regime V ^ 0. This is illustrated in Fig. [4]d. 
The high efficiency regions lie within the dot resonances 
(indicated by the dashed white lines) which are shifted 
from zero energy due to the presence of an induced mini- 
gap by the proximity with the superconducting lead. 

The origin of these high efficiency regions can be un- 
derstood qualitatively from the spectral density on each 
dot, which is shown in Fig. [dj;. At any given point along 



the line El ~ —En electron and hole states are split 
due to hybridization between the dots. Furthermore, the 
electron-hole symmetry in the local spectral density is 
lost along this line. Crossed Andreev processes like the 
one sketched as the CAR arrow in Fig. |4j; combines elec- 
tron and hole states on each dot with high spectral den- 
sity. These inter-dot transitions are then more favorable 
than the intra-dot electron-hole conversions (arrows AR 
and AL), in which either the electron or the hole state 
has low spectral density. As a consequence, local An- 
dreev processes become suppressed while non-local CAR 
processes are enhanced, thus explaining the efficiency in- 
crease. 

Conclusions: We have analyzed the splitting efficiency 
of SWCNT double quantum dot devices in terms of ma- 
terial and geometrical parameters. The single channel 
character of the connection between the dots in this con- 
figuration explains the weak decay of CAR with distance 
which is consistent with the experimental observations. 
Furthermore we have shown how the splitting efficiency 
can rise up to ~ 100% by working in the non-linear 
regime. We expect that our analysis can guide future ex- 
periments for the production of entangled electron pairs 
using these devices . 
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